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A Floquet scattering approach to parametric electron pumps is presented and compared with 
Brouwer's adiabatic scattering approach [Phys. Rev. B 58, R10135 (1998)] for a simple scattering 
model with two harmonically oscillating (5-function barriers. For small strength of oscillating po- 
6J[), tentials these two approaches give exactly equivalent results while for large strength, these clearly 

deviate from each other. The validity of the adiabatic theory is also discussed by using the Wigner 
delay time obtained from the Floquet scattering matrix. 

• PACS number(s): 73.23.-b, 72.10.Bg, 73.50.Pz 

J* ; I. INTRODUCTION 

A parametric electron pump is a device that generates a dc current at zero bias potential through cyclic change 
of system parameters. The most direct way to create a dc current was originally proposed by Thouless Q], who 
considered a system subjected to a traveling wave potential. This can be realized for example with the help of 
surface acoustic waves [^). Another possibility is to utilize quantum dots. In closed systems operating in the Coulomb 
blockade regime, integer number of electric charge can be transferred by sequential changes of barriers like a turnstile 
H , whereas in open systems the electron pumping can be driven by adiabatic shape change in the confining potential 
or other parameters which affect the interference pattern of the coherent electrons in the device. After a cycle of the 
adiabatic shape change we return to the initial configuration, but the wavefunction may have its phase changed from 
O , the initial wavefunction. This is the so-called geometric or Berry's phase B. This additional phase is equivalent to 
some charges that pass through the quantum dot, namely, pumped charge ^|,^|. Recently, adiabatic charge pumping 
■r— ( | in open quantum dots has attracted considerable attention [p[-p^| . 

Switkes et al. report an experimental investigation of electron pumping through an open quantum dot under the 
shape deformation controlled by two ac gate voltages For weak pumping the dc voltage induced by the pumped 
dc current not only have a sinusoidal dependence on phase difference <f> of the two ac voltages applied to the gates, 
qq , but is also proportional to the square of the pumping strength A. Many aspects of the experimental results can be 
understood in terms of Brouwer's scattering approach JtJ. The change in the charge of the dot in response to a small 
s ! ! variation of external parameter SX is given by SQ = SXidn/Xi, where dn/dX is the emissivity fl3| . The pumped 
charge during each cycle can be determined by integrating SQ along the closed path in the parameter space {X;} 
defined by one pumping cycle. Then, the pumped charge Q is rewritten as an integral over the surface enclosed by 
the path by using Green's theorem, e.g. in two parameter case Q = e J A dX 1 dA" 2 n(AT 1 , A" 2 ). Assuming II is constant 
over small area of (Xi,X%), for weak pumping one can show the pumping current I p oc wA 2 sin0, which exactly 
corresponds to what was observed in the experiment. For strong pumping, however, the dc voltage generated by 
the pumped current deviates from the sinusoidal dependence on 0, and also departs from the A 2 dependence. The 
non-sinusoidal dependence can be understood by taking into account that II is no longer constant over the integration 
area when A is considerably large. In Ref. ||] the anomalous A dependence was ascribed to the occurrence of significant 
heating and dephasing as a result of strong pumping. 

The parametric electron pump is a time dependent system driven by (at least) two different time periodic per- 
turbations with the same angular frequency and phase difference (f>. One can deal with this problem using not only 
adiabatic approximation exploited by Brouwer but also Floquet approach Jl4| . An oscillating potential can transfer an 
incoming electron of energy E to Floquet side bands at E ± nfiu, where n is an integer and uj is the angular frequency 
of the oscillation. A scattering matrix for a time dependent system can be constructed from the interplay of these 
sidebands p5[ , p!6[ . So far there have been relatively few works on non-adiabatic parametric electron pumps p7| , ^8| , 
and the comparison between the adiabatic and the Floquet approach is still missing. We show, using a simple scat- 
tering model with two harmonically oscillating ^-function barriers, that for small A these two approaches give exactly 
equivalent results for the pumped charge while for large A they are different. Even though the pumped current as a 
function of A obtained from these two approaches shows qualitatively similar behavior, the physical interpretation is 
completely different. This discrepancy in large A is not fully ascribed to the breakdown of adiabatic conditions since 
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the Wigner delay times calculated by using Floquet scattering matrix are still much smaller than an oscillation period 
of external pumping potentials. 

In Sec. II, we introduce our model system and investigate the characteristics of the pumped current using Brouwer's 
scattering approach. In Sec. Ill, we study the Floquet approach for parametric electron pumps and compare with 
those of Sec. II. We also discuss the validity of adiabatic theory by using the Wigner delay time. Finally, we conclude 
our paper in Sec. IV. 



II. BROUWER'S SCATTERING MATRIX APPROACH 



As a model system we consider ID two harmonically oscillating (^-function barriers with the strengthes X\ = 
V\ + Ai cos lot and X 2 = V% + A2 cos(wt + </>) respectively, separated by a distance d. This is a simplified model of 
the experiment by Switkes et al., but possesses many important characteristics and can be easily handled. Wei et 
al. studied parametric charge pumping aided by quantum resonance using this model and found that the pumped 
current has large values near a resonance level [|l9| . Due to the double barrier geometry, resonant tunneling also plays 
an important role in charge pumping. 

The 2x2 scattering matrix of the double barriers with the strengthes X\ and X2 is given by 

"-(It). 0) 



J r' / 

where r and t are the reflection and the transmission amplitudes respectively, for modes incident from the left; r' and 
t' are similar quantities for modes incident from the right. The charge emitted per cycle to the right is obtained from 







dni dX\ dni dX 2 



where 



dn ^ :m ( JL t * + ^r>* ] (3) 



dX m \ dX m dX n 

(m = 1,2), and T p (— 2tt/u)) is the period of the pumping. One can show that the emitted charge to the left Q 2 is 
equal to —Q\. Equation (0) can be rewritten in the following form by using Green's theorem 



Qi = e 



[ dx 1 dx 2 n(x 1 ,x 2 ), (4) 

J A 



where II(Xi, X2) = d(dn\ / dX 2 ) / dX\ — d{dn\/ dX\) / dX 2 . The pumped current is easily obtained from I\ = Q\/T p . 

We use the parameters d = 50 nm for the distance between the two barriers, the effective mass /i = 0.067m e of an 
electron in GaAs, and T p = 9.09 ps for the period of pumping, which corresponds to Tilo — 0.45 meV. Figure [l] shows 
the pumped current as a function of the energy of an incident electron with V% — V 2 — V , Ai = A 2 = A, and <f> = tt/2. 
We present two examples; the first is a nearly open case (V=0) in Figs. [l]a and [ljj, and the second is a closed case 
(V—225 meV-nm) in Figs, ^jp and |l|d. All of them show interesting resonance-like structures. In the closed case one 
can directly see the relation between the pumped currents and transmission resonances in Figs, ^c and [l]d, where the 
transmission poles are denoted by the filled circles (•). To find these poles we considered a scattering problem with 
static double barriers of strength V. The imaginary part of the pole (denoted by the length of error bars attached 
to the filled circles) is related to the lifetime of a resonant (or quasibound) state, and determines the width of the 
transmission peaks. Peaks of the pumped current and their width are also described by the real and imaginary part of 
the resonance poles respectively, as shown in Fig. [I]c. In Fig. ^d, however, the width of the pumped current is larger 
than the imaginary part of the resonance energy due to the effect of strong oscillation. Moreover, in the open case 
(Figs. 0a and |l|b), this analysis is no more relevant since there is no resonance in the static double barrier model. 

Figure ^| shows the pumped current as a function of the pumping strength A for V = (the open case) with cf> = tt/2. 
Here for small A the pumped current depends on A 2 and has an exact sinusoidal dependence on <fi, while for larger A 
the pumped current saturates and even decreases as A increases, and also deviates from the sinusoidal dependence on 
4>. The decrease of the pumped current can be understood by considering the charge flux H(Xi, X 2 ) in Fig. |^. For 
small A, the integration area only contains the positive IPs (the solid circle). Hence, as A increases J. II also becomes 
larger, which is roughly proportional to the area enclosed by the circle, i.e. J A II oc A 2 . Around A = 225 meV-nm 
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(the dashed circle in Fig. ||) the integration area begins to include the region with both the positive and the negative 
ITs, which means J A II no longer increases as A increases. This explains why the pumped current saturates and even 
decreases above A = 225 meV-nm as shown in Fig. ||. From this we expect that the behavior of the saturation of the 
pumped current would depend on (ft since the shape of the integration area is determined by ip. For example, with 
<f> = 1.1-k the integration area denoted by the dashed dotted curve in Fig. || is distorted from the original circular 
shape, so that even for A = 675 meV-nm the integration area contains only positive ITs. Figure || clearly shows the 
pumped current with (f> = l.ln saturates much more slowly in comparison with the pumped current with <f> = Tt/2. 
The non-sinusoidal dependence of the pumped current on (f> for large A is also ascribed to the loss of circular symmetry 
of II in (X 1 ,X 2 ). 



III. FLOQUET APPROACH 



Now we study the Floquet approach of the problem investigated in Sec. II. Using the scattering matrix of a single 
5- function with sinusoidal time dependence (see Appendix for details), we can obtain the total scattering matrix of 
the oscillating double ^-functions in the following form 

^=(w)' (5) 

where r = ri+t^(I — Qr i? Qr L ) _1 Qr i? Qt i , and t = t^(I — Qr L Qr fl ) _1 Qt £ ; r' and t' can also be obtained by 
replacing L by R in r and t respectively. Here ^r(l) an d ^r(l) are the reflection and the transmission matrices 
respectively, for the right (left) delta function with time dependence for modes from the left, and I is an identity 
matrix. Due to the reflection symmetry of each delta function t'r^ = ?r(l), and ^'r(l) = ^R(L)- During each one-way 
trip an electron at energy E n — E + nhu picks up a phase factor exp(ik n d) , which is represented by the diagonal 
matrix Q„ m = enp(ik m d)S mn - From Sfi we can obtain the total transmission coefficients for the propagating mode 
entering in the mth channel 

oo 

T^ m {Epi) = 2_j ; (6) 

n=0 

where Epi is the Floquet energy. The total transmission coefficient from the left to the right as a function of an 
energy of an incident electron is given by T^(E) = T^ m (Epi) where E = Epi + mhuj. The total transmission from 
the right to the left T^(E) can also be determined in a similar way. 



The pumped current to the right is given by [[17 20 



h^f dEdE'[t(E\ E)f L (E) - t'(E, E')f R (E% (7) 

where t(E',E) represents the transmission probability for scattering states incident from the left at energy E and 
emerging to the right at E' , and t'(E,E') is defined in a similar manner for the reverse direction, /x (Jr) is the 
Fermi-Dirac distribution in the left (right) reservoir. Since the Floquet energy must be conserved during the scattering 
process, E' and E are given by Epi + nhu and Epi + rriTiuj respectively (n and m are integers). Using Eq. (|^), we 
can derive the relation J dE't(E', E) = J2 n ^nm(^Fi) — T^{E). Without external bias (/l = Sr = f) Eq. (0) can be 
rewritten as 



h = \\ dEf[E)[T . ( E ) T LEi], (S) 



At zero temperature it becomes 



2e 



h = -j dE[T^(E)-T^(E% (9) 



where Ep is the Fermi energy. Equation (g) can also be expressed as I\ =/_,. — J<_. This is quite interesting because 
from this point of view the pumped current merely corresponds to the difference of two currents having the opposite 
directions through a scatterer. 
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Figure [I] shows that in the open case (V = 0) for small A (Fig. pi) the pumped current obtained from the Floquet 
approach is equivalent to that of Brouwer's while for rather larger A (Fig. |l|b) they deviate from each other. When 
V 7^ (the closed case), even for small A (Fig. |]c) they are quantitatively different near the resonances. Figure || 
shows TU and T*_ in comparison with the pumped current I±. In the closed case (Figs. [|c and ^d) it is shown that 
TU ~ T<— , and the transmission resonances are directly related to the maxima of the pumped currents. In the open 
case (Figs. ||a and ^b), however, and TU considerably differ from each other, and it is hard to determine the 
relation between the resonances and the pumped currents. In fact, looks out of phase with T<_, and their resonance 



structures are quite complicated. In the open case the resonance-assisted electron pumping |19 seems unclear even 
though the oscillatory behavior of transmission looks similar to that of the pumped current. 

As A increases with V = 0, in Fig. |^, the pumped current obtained from the Floquet approach also saturates and 
even decreases although their exact values are different from those of Brouwer's. We also plot 7_> and /<_ in Fig. |^, 
where it is clear that in the Floquet approach the decrease of I-y for large A is ascribed to the decrease of both 7_^ and 
7<_. Usually the stronger the barrier strength, the smaller the transmission. This interpretation differs from that of 
Brouwer's, in which the decrease of the pumped current is explained by considering the structure of IT(Xl,X2) (Sec. 
II). Since we cannot define H{X\,X2) explicitly in the Floquet formalism, we plot II integrated over small circular 
area with a radius A centered at (Vi, Va). We consider only the case that V\ = V% = V. Figure [| shows that the 
integrated IPs obtained from the two different approaches for A = 22.5 meV-nm look very similar to each other, and 
their overall structure are governed by the static double barrier resonances denoted by the circles (o). 

Brouwer's approach is based upon the adiabatic approximation, which implies that any time scale of the problem 
considered, especially the electron dwell time in a quantum dot (or inside double barriers in our case), must be much 
smaller than the period of the oscillation of a external pumping T p |?J . Using the Floquet formalism we can calculate 
the Wigner delay time tw , which is the interaction time of the incident electron with the scattering potential pl[ (see 
also p2| , p3| |). In this sense tw corresponds to the electron dwell time in the quantum pump. To obtain the Wigner 
delay time we use the eigenvalues of the scattering matrix Sfi- Due to the unitarity of Sfi all the eigenvalues lie on 
the unit circle and can be written in the form exp(i8 a ). The Wigner delay time is defined by 

Tw = hY,^§\{k n \e a )\\ (10) 

a 

where the eigenstate corresponding to the eigenvalue 9 a and an input propagating state (or channel) with momentum 
k n are denoted by \9 a ) and \k n ) respectively p4| . It is worth noting that the Wigner delay time is a function of 
the energy of the incident particle E (= Epi + nfiu)), and \k n ) and \6 a ) are determined by n and Epi respectively. 



If (k n \6 a ) is ignored in Eq. (10) the Wigner delay time tw becomes trivial, i.e. tw{E + nfiuj) — tw(E). Thus, we 
cannot observe any signature of the double barrier resonances. 

Figure |^ shows the Wigner delay time using the same parameters exploited in Fig. |5|. The Wigner delay times 
become smaller for both larger E and lower V, which can be understood when we take into account that usually the 
electron dwell time is short if the energy of an incident electron is large or the scattering barrier is weak. The Wigner 
delay times have larger values near the resonances, which is ascribed to the fact that at the resonances an electron can 
stay in the quantum dot for a long time. Near the resonances the adiabatic condition can break down. This explains 
the deviation between the pumped currents observed in Figs. p. In the open case (V = 0) we also check the Wigner 
delay time for rather larger A (up to 675 meV-nm), and observe they are smaller than T p by two orders of magnitude 
except for small incident energy E (not shown), which means Brouwer's approach should be still applicable even for 
large A. This leads us to conclude that the deviations of the pumped current observed in Figs. 0a, [ij) and || for large 
A are not simply ascribed to the breakdown of adiabatic condition. It is worth noting that recently it is shown that 
Brouwer's formula is valid for arbitrary amplitudes of the modulating potential, as long as the lowest-order adiabatic 
approximation can be employed [ jlj) . Let us just note that in Fig. ^ the discrepancy for large A appears when |7i |/|7_^| 
(or |7i|/|7,_|) is not so small. 

One of the interesting consequences from Eq. (Q) is that the pumped current still exists even in the cases cf> — or 
7r when Ai 7^ A2 fTs) . Since the integration area in parameter space is zero when cf> = or it, in Brouwer's approach 
the pumped current definitely vanishes. In contrast, even when <f> — or 7r, an asymmetry of the potential can lead to 
the asymmetry of the currents |2(| , which is nothing but the pumped current in Eq. (||) . Figure shows the pumed 
current as a function of ratio of the strength of two barriers A2/A1 with <p = or ir, and E = 6.005?ia;. Note that the 
pumped current is zero when Ai = A2. The oscillatory behavior is also related to the double barrier resonances. 
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IV. SUMMARY 



We investigate the Floquet scattering in parametric electron pumps in comparison with Brouwer's adiabatic scat- 
tering approach exploiting two harmonically oscillating S- function barriers. In the Floquet approach the pumped 
current simply corresponds to the difference of two currents having the opposite direction through a scatterer. For 
small strength of the oscillating potentials, these two kinds of approach give exactly equivalent results while for large 
strength these show deviation. Even though for large strength we obtain qualitatively similar results for the pumped 
currents using both approaches, the physical interpretation is completely different. The validity of the adiabatic 
approximation is also discussed by calculating Wigner delay time. For large A although the adiabatic condition is well 
satisfied, a quantitative discrepancy between both approaches is still observed. In the Floquet approach, a non-zero 
pumped current can be obtained even when <j> — or 7r (no current at all in Brouwer's approach), if the spatial 
symmetry of the potential is broken (Ai ^ A2). 
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APPENDIX: 



The scattering problem of a single 5-function impurity with sinusoidal time dependence has been investigated by 
several authors [p5|-p7f . We would like to summarize how to construct its Floquet scattering matrix in this Appendix. 
The system is described by the Hamiltonian 



H(x,t) 



2/i dx 2 



+ [V s + V d cos(ujt + <j))]6(x), 



(Al) 



where fi is the mass of the incident particle, while V s and Vd represent the strength of the static and the oscillating 
potential respectively. Using the Floquet formalism the solution of this Hamiltonian can be expressed as 



(A2) 



where Epi is the Floquet energy which take continuous values in the interval < Epi < Tuj. 
Since the potential is zero everywhere except at x = 0, ip n {x) is given by the following form 



A n e lk ™ x + B n e- iknX , x<0 

Q ik n X + D e -ik„X X>Q 



where k n — y / 2/i(Epi + nfiui) /Ti. The wave function ^e fi (x, t) is continuous at x = 0, 

A n + B n = C n + D„, 

and the derivative jumps by 



dx 



Ep-i 



x=0+ 



dx 



x=Q- 



2m 



2 [V 3 + V d CO S (LOt + 4>)]^E Fl (0, t). 



Using Eq. (A2) this leads to the condition 



(A3) 



(A4) 



(A5) 



ik n (C n — D n — A n + B n ) 
= ls {A n + B n ) + ld (e-^A n+1 + e^K-i + e~^B n+1 + e i ^ n _ 1 ) (A6) 
= 7s (C„ + D n ) + 7 d(e- # C n+1 + e**C n _i + e^D n+1 + e*i? n _i), 

where 7 S = 2fj,V s /h 2 and 7^ = fiVd/h 2 . After some algebra we have the following equation from Eqs. (A4) and (A6) 
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where 



-(j + r)-^ 
(i + r)- 1 



(i + r)- 1 



(A7) 



V 










7 s /ifc_i 








Ide^/iki 







'Js/ikl Jd£~ 



b /ih 



(A8) 



and / is an infinite-dimensional square identity matrix. Eq. (A7) can also be expressed in the form |out) = M |in), 
where M connects the input coefficients to the output coefficients including the associated evanescent Floquet side- 
bands. In order to construct the scattering matrix we multiply an identity to both sides, K~ X K |out) = MK^ 1 K |in), 
where K nm = \fk~^8 nm . Then we have J ou t — AlJi n , where J represents the amplitude of probability flux and 
M = KMK -1 . It should be mentioned that M is not unitary due to the evanescent modes included. 

If we keep only the propagating modes, we obtain the unitary scattering matrix S Jj5| , fL6| |, which can be expressed 
in the following form 



/ roo r i 
rio rn 



S : 




c 00 


r 


• • • f' 


t'n •• 


■ ' ' r 00 


fit ■ 







(A9) 



where r nm and t nm are the reflection and the transmission amplitudes respectively, for modes incident from the left; 
r' nm and t' nm are similar quantities for modes incident from the right. 
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FIG. 1. The pumped current 7i calculated by using Brouwer's approach (solid curves) and the Floquet approach (dotted 
curves) with <f) = n/2 for (a) A = 22.5 meV-nm and (b) A = 225 meV-nm with V = 0, and (c) A = 22.5 meV-nm and (d) A = 225 
meV-nm with V = 225 meV-nm. In (c) and (d) the transmission resonances are denoted by the filled circles (•) obtained from 
considering static double barriers, whose y values are chosen arbitrarily. The attached error bars represent the sizes of the 
imaginary energy of each resonance. 
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FIG. 2. The pumped current Ji as a function of A with Vi = V2 = and i5 = 6.005?k<; by using Brouwer's approach for 
<j> — 7r/2 (the solid curve) and (f> — l.ln (the dashed curve), and the Floquet approach for 4> = tv/2 (•) and cf> — l.ln (o). 
x and + represent and 7^ for <f> — tt/2 respectively. The inset shows the pumped current h as a function of <f> for (a) 
A = 22.5 meV-nm and (b) A = 675 meV-nm by using Brouwer's approach (the solid curves) and the Floquet approach (the 
dotted curves). 
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FIG. 3. Contour plot of flux H(Xi, X2) for E = 6.005?icj. The solid and the dashed contours represent positive and negative 
fluxes respectively, in the range of —0.016 to 0.016. The circles represent the integration area for <j) = n/2 (see the text for 
detail) with A = 22.5 meV-nm (thick solid curve), A = 225 meV-nm (thick dashed curve), and A = 675 meV-nm (thick dotted 
curve). The dashed dotted curve represents the integration area with <f> = 1.1-7T and A = 675 meV-nm. 
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FIG. 4. Pumped current (thick solid curves), TU (thin solid curves), and T«_ (dotted curves) with the same parameters as 
used in Fig. [[]. For clear comparison the unit and the scale of the y-axis are arbitrarily chosen. 
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FIG. 5. Contour plot of the pumped current Ii obtained from integrating II with A = 22.5 meV-nm (see the text for detail) by 
using (a) Brouwer's approach and (b) the Floquet approach. The solid and the dashed contours represent positive and negative 
currents respectively, in the range of —0.2 nA to 0.2 nA. The open circles (o) and the vertical lines represent the resonance 
energies of quasibound states of static double 5-function barriers and the resonance energies from the condition ir Ti m /2(id 
with an integer m (m = 1, 2, ■ ■ •) respectively. 
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FIG. 6. Grey scale plot of the Wigner delay time. Black denotes times larger than 0.1T P . The open circles (o) represent the 
resonance energies which are the same as in Fig. [| 
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FIG. 7. Pumped current I\ as a function of potential asymmetry A2/A1 with Ai = 22.5 meV-nm and E = 6.0057i(jj at (f> = 
(solid curve) and <j> = n (dashed curve), where the currents are normalized to their values at at = tt/2. Inset shows the 
magnification of a part of the plot. 
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